Correlation Length Exponent in the Three-Dimensional Fuse Network 
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We present numerical measurements of the critical correlation length exponent v in the three- 
dimensional fuse model. Using sufficiently broad threshold distributions to ensure that the system is 
the strong-disorder regime, we determine v tohe v — 0.86 ±0.06 based on analyzing the fluctuations 
of the survival probability. The value we find for f is very close to the percolation value 0.88 and we 
propose that the three-dimensional fuse model is in the universality class of ordinary percolation. 



It is already twenty years since the publication of the 
first experimental evidence of scaling in the morphology 
of brittle fractures 0. About seven years later it was 
proposed that not only is there scaling, but the scaling 
properties are universal, in the sense that they do not de- 
pend on material properties 0, Q • There is now mount- 
ing evidence for this hypothesis, which may be expressed 
as the scaling invariance 



t:{z; X, y) — A^7r(A^z; Aa;, Xy) 



(1) 



of TT, which is the probability density that at position 
{x, y) in the average fracture plane, the fracture is at 
height z given that it is at z = at (0, 0), with ( as the 
universal roughness exponent having a value very close 
to 0.80 for a large class of materials. One experimentally 
important consequence of this scaling is that the average 
fracture width w scales as 



w 



(2) 



where L is the linear size of the average fracture plane. 

Ever since the proposal of universality, it has remained 
a theoretical challenge to explain this value. Recently, it 
was suggested by Hansen and Schmittbuhl that it has its 
origin in the fracture process being a a correlated perco- 
lation process 0|. The essence of the argument is based 
on existence of a localization length I and a correlation 
length ^ that grows during the breakdown process. The 
localization length depends on the disorder in the mate- 
rial: Stronger disorder means larger localization length. 
Whether the localization length diverges for large but 
finite disorder or it only reaches this limit for infinite dis- 
order is at present not known. However, mean field argu- 
ments suggest that the former scenario is the correct one 
For correlation lengths ^ much smaller than the lo- 
calization length I, Hansen and Schmittbuhl Q assumed 
a relation 



(3) 



where p is the local damage density and pc is the damage 
density at failure. This relation is taken directly from 
percolation theory. The reason it is only valid for large 
localization lengths I is that p is assumed to be spatially 
stationary (meaning that the statistical distribution of p- 
values is independent of position). The correlation length 
exponent v has the value 4/3 in two-dimensionalperco- 
lation and 0.88 in three-dimensional percolation It is 
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Figure 1: Normalized histogram where each bin is averaged 
over a plane orthogonal to the (l,l,l)-direction for D = 10. 



by no means given that v should be the same in the brit- 
tle fracture problem — and Toussaint and Pride suggest 
that it is equal to 2 Q. However, it was suggested by 
Hansen and Schmittbuhl that the two-dimensional fuse 
model has = 4/3 placing it in the same universality 
class as two-dimensional percolation. When the correla- 
tion length approaches the localization length /, gradients 
develop in the damage — p can no longer be regarded as 
spatially stationary — and using arguments from gradi- 
ent percolation Hansen and Schmittbuhl suggested 
the relation 



c 



2v 



l + 2iy 



(4) 



With z/ = 4/3 for the two-dimensional fuse model, this 
leads to C = 8/11 « 0.73. Recent numerical calculations 
gives C = 0.74 ±0.03 Q. 

Recently, Kumar et al. ^3 have proposed that there 
is no universal correlation length exponent v in the two- 
dimensional fuse network. The numerical evidence pre- 
sented is based on a disorder having a small, finite local- 
ization length so that p is not spatially stationary due 
to localization. However, the analysis implicitly assumes 
that Eq. JSJ is valid, which requires p to be spatially sta- 
tionary. Hence, there is no support for the conclusion 
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Figure 2: Log-log plot of the fluctuations of the density of 
broken bonds Wc = ((p^) — (p)^y^^ against L. The disorder 
is hence D = 10, 12, 15 and 20 respectively. The slopes are 
for D = 10: 1.11, D = 12: 1.14, D = 15: 1.10 and D = 20: 
1.23. Their mean is 1/v = 1.16±0.06, giving v = 0.86±0.06. 



reached. 

It is the aim of this letter to measure v in the three- 
dimensional fuse model. We find the value v = 0.86 ± 
0.06. This is close to the three-dimensional percolation 
value V — 0.88, hence supporting the notion that the fuse 
model is in the universality class of ordinary percolation, 
both in two and three dimensions. The roughness ex- 
ponent C was measured by Batrouni and Hansen to 
be C = 0.62 ± 0.05. Using Eq. gj) with v = 0.86, we 
find C = 0.63. Hence, the value for v we report here 
is consistent with the roughness exponent measured in 
when using Eq. We note, however, that this 

value for the roughness exponent is not consistent with 
the one reported by Raisanen et al. who reported a 
roughness exponent close to the minimal energy result, 
C = 0.41 ± 0.02 T^l, claiming that they should be identi- 
cal. 

The fuse model that we study consists of an oriented 
simple cubic lattice. As in Ref. pTl | , w e use periodic 
boundary conditions in all directions |l4| and the aver- 
age current flows in the (l,l,l)-direction. Each bond is an 
ohmic resistor up to a threshold value. When this value is 
reached, the resistor turns irreversibly into an insulator. 
The threshold values are drawn from a spatially uncor- 
related probability density p{t). A voltage drop equal to 
unity is set up across the lattice along a given plane or- 
thogonal to the (l,l,l)-direction. The currents are then 
calculated using the Conjugate Gradient algorithm [T5l| . 
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Figure 3: (p) plotted against L~^'" with v = 0.88. D = 
10 (x), D = 12 (A), i3 = 15 (*) and D = 20 (+). As L -> oo 
the straight lines extrapolate to the thresholds pc — 0.573, 
Pc = 0.587, Pc = 0.602 and pc = 0.616 respectively. 



After the currents i have been determined, the bond hav- 
ing the largest ratio max(i)/max(t) is determined. This 
bond is then removed and the currents are recalculated. 
We do not allow the final crack to cross the plane along 
which the voltage drop is imposed. This simplifies the 
analysis of the final crack breaking the network apart, 
while it only imposes weak finite size corrections to frac- 
ture patterns. 

The threshold values t constructed by setting t ~ r^ , 
where r is drawn from a uniform distribution on the unit 
interval |l6j| . This corresponds to a probability density 
p{t) cx t-^+^ on the interval < i < 1 with P = 1/D. 
The parameter D > controls the width of the distribu- 
tion: Larger values of D corresponds to stronger disor- 
der. In order to ensure that our results are obtained in 
the strong disorder phase of the fuse model, we studied 
D — 10, 12, 15 and 20. Our system sizes varied from 
i = 6 to 24 with 5000 samples generated for the smallest 
sizes to 200 samples for the largest sizes. 

With D — 20, the smallest threshold values gener- 
ated are of the order (24^)"^'^ w 10"^^. The system 
has, however, still not entered purely screened percola- 
tion regime. With this level of disorder, the system fails 
when a fraction of about 0.62 of the bonds have failed. 
The threshold values of the bonds that fail near the end 
of the process are about 0.62^" « lO""* — which is of the 
order of the currents that are carried by the bonds in the 
system. Hence, there is competition between threshold 
values and currents, making the failure process a corre- 
lated one rather than a pure percolation one even in this 
seemingly extreme case. 

Fig. n shows the damage profile in the current direc- 
tion of the random fuse model with D — 10. We de- 
note the (l,l,l)-direction the z-direction. We define the 
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Figure 4: pc plotted against 1/D and extrapolated to infi- 
nite disorder giving Pc(oo) = 0.66. Extrapoating the straight 
line, towards increasing l/D- values, we find that pc becomes 
negative for 1/D > 0.75. 

damage as the normalized average number of burned-out 
fuses in the plane orthogonal to the z-direction at z. The 
distribution has a weak maximum in the middle. This 
indicates a finite but large localization length Such a 
maximum is smaller or entirely absent from the stronger 
disorders (i.e. larger Z)-values) we studied. 

Following percolation analysis 0, we define the sur- 
vival probability 11 indicating the relative number of lat- 
tices that has survived for a given average damage p. 
Assuming that the disorder is broad enough so that p 
is independent of z and there is a finite critical value of 
p = Pc at which 50 % of the lattices survives, we have 
that 

n = n[(p-p,)Li/^]. (5) 

This scaling ansatz implies that both the mean value 
of the density of broken bonds (p) and the fluctuations 
((P^) ~ {p)^y^^ at breakdown scales as L"'^/'^ using 




and 

In Fig. 121 the fluctuation s of the density of broken 
bonds, Wc = \/ {p^) — (p)^, have been plotted against 
the system size L. The mean value of the slopes gives 



V = 0.86 ±0.06 which is consistent with percolation value 
1/ = 0.88. 

Using v = 0.88 from standard percolation we now turn 
to the scaling of {p). From finite-size scaling analysis, we 
expect the functional dependency 

{p)=Pc-^ (8) 

on L. We show this relation for different values of D in 
Fig. El 

This way of measuring the critical exponent v is much 
less sensitive than the one presented in Fig. El From 
standard percolation in a simple cubic lattice the thresh- 
old for an infinite system is Pc = 0.752 '6^. The ex- 
trapolations done in Fig. ^ show results lying below this 
threshold. However, this is to be expected as the per- 
colation process in this limit is screened |17|| . This re- 
sult strongly indicates that there is a strong disorder 
regime for finite disorders with pc larger than zero in 
the three-dimensional fuse model. In fact, extrapolat- 
ing the straight line in Fig. ^ towards larger 1 / D- values 
will result in pc reaching zero and becoming negative at 
D < 1.33. This is physically impossible and Pc remains 
zero in this range. This indicates that there is a transition 
from a percolation-like regime with pc > for D > 1.33 
to a regime with Pc ^ for D < 1.33. This latter regime 
has been described as the diffuse localization regime in 
i. 

In summary, we have determined the correlation length 
exponent in the three-dimensional fuse model to be = 
0.86 ± 0.06. This is consistent with the percolation value 
oi — 0.88. Furthermore, using Eq. Q), this is consis- 
tent with the previously measured roughness exponent 
C = 0.62 ±0.05 lending support to the scenario pro- 
posed by Hansen and Schmittbuhl 4] for understanding 
the universality of the roughness exponent in the fuse 
model and brittle fracture. Our analysis was based on 
studying the fuse model with strong enough disorder for 
the breakdown process to develop in a percolation-like 
manner with p spationally stationary so that the tools 
developed for studying that problem could be used in the 
present one. We note that in this regime, one will not 
see the fracture roughness scaling of Eq. (Q): The frac- 
ture will have a fractal structure. When, on the other 
hand, the disorder is weak enough for localization to set 
in, p is no longer spatially stationary, making a direct 
measurement of v based on fluctuations in p impossible. 
However, it is in this regime fracture roughness scaling 
as in Eq. is seen as shown in .IjJ. 
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